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Abstract 

The spin structure functions g\ and gi have been calculated in the resonance 
region and for small and intermediate momentum transfer. The calculation 
is based on a gauge-invariant and unitary model for one-pion photo- and 
electroproduction. The predictions of the model agree with the asymmetries 
and the spin sturcture functions recently measured at SLAC, and the first 
moments of the calculated spin structure functions fullfil the Gerasimov-Drell- 
Hearn and Burkhardt-Cottingham sum rules within an error of typically 5- 
10 %. 
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I. INTRODUCTION 



The spin structure of the nucleon has been studied by increasingly accurate measurements 
since the end of the 70's. By scattering polarized lepton beams off polarized targets, it has 
become possible to determine the spin structure functions gi(x, Q 2 ) and g2(x,Q 2 ). These 
functions depend on the fractional momentum of the proton, x = Q 2 /2mz/, where — Q 2 is the 
square of four-momentum transfer, v the energy transfer in the lab frame, and m the nucleon 
mass. The results of the first experiments at CERN [|IJ and SLAC sparked considerable 
interest in the community, because the first moment of gi, V\ = Jq g\{x)dx, was found to be 
substantially smaller than expected from the quark model, in particular from the Ellis- Jaffe 
sum rule || . This "spin crisis" of the nucleon led to the conclusion that only a small fraction 
of the nucleon's spin is carried by the valence quarks, and that the main spin contributions 
should be due to sea quarks, quark orbital momentum and gluons. However, the difference 
of the proton and neutron moments, T'l — T", was found to be well described by Bjorken's 
sum rule |4]], which is a strict QCD prediction. 

Recent improvements in polarized beam and target techniques have made it possible 
to determine both spin structure functions over an increased range of kinematical values. 
The reader is referred to the publications of the collaborations at CERN ||, SLAC || 
and DESY J7| and to the references on earlier work quoted therein. While most of these 
investigations have concentrated on the high Q 2 region, the recent measurements of the 
E143 Collaboration at SLAC f| have been extended to the range of 0.6 (GeV/c) 2 < Q 2 < 
1.2 (GeV/c) 2 , which bridges the gap between low-Q 2 nonperturbative resonance phenomena 
and high-Q 2 perturbative QCD. 

It is the aim of our contribution to present the results of the recently developed Unitary 
Isobar Model (UIM, Ref. ||) for the spin asymmetries, structure functions and relevant 
sum rules in the resonance region. This model describes the presently available data for 
single-pion photo- and electroproduction up to a total cm energy W mSuX ~ 2 GeV and for 
Q 2 < 2 (GeV/c) 2 . It is based on effective Lagrangians for Born terms (background) and 



2 



resonance contributions, and the respective multipoles are constructed in a gauge-invariant 



and unitary way for each partial wave. The eta production is included in a similar way [JTO 
while the contribution of more-pion and higher channels is modeled by comparison with the 
total cross sections and simple phenomenological assumptions. 

The spin structure in the resonance region has intrigued many authors because of an 
expected rapid change from the physics of resonance-driven coherent processes to incoherent 
scattering off the constituents. In particular the first moment Ti is constrained, in the limit 
of Q 2 — > (real photons), by the famous Gerasimov-Drell-Hearn sum rule (GDH, Ref. [|TT|] ), 
Ti - * —Q 2 K 2 /8m 2 , where k is the anomalous magnetic moment of the nucleon. The reader 
should note that here and in the following we have included the inelastic contribution to 1^ 



only. As has been pointed out by Ji and Melnitchouk ||12|| , the elastic contribution is in fact 
the dominant one at small Q 2 and has to be taken into account in comparing with twist 
expansions about the deep- inelastic limit at those values of Q 2 . 

In the case of the proton, the GDH sum rule predicts 1^ < for small Q 2 , while all 
experiments for Q 2 > 1 (GeV/c) 2 yield positive values. Clearly, the value of the sum rule has 
to change rapidly at low Q 2 , with some zero-crossing at Q 2 , < 1 (GeV/c) 2 . The evolution of 
the sum rule was first described by Anselmino et al. [HJ in terms of a parametrization based 
on vector meson dominance. Burkert, Ioffe and others Hl4| , |T5| refined this model considerably 
by treating the contributions of the resonances explicitly. Soffer and Teryaev [|]1J suggested 
that the rapid fluctuation of Ti should be analyzed in conjunction with T 2 , the first moment 
of the second (transverse) spin function. The latter is constrained by the less familiar 
Burkhardt-Cottingham sum rule (BC, Ref. JT7j) at all values of Q 2 . Therefore the sum of 



the two moments, Ti + I^, is known for both Q 2 = and Q 2 — > oo. Though this sum 
is related to the practically unknown longitudinal-transverse interference cross section o' LT 
and therefore not yet determined directly, it can be extrapolated smoothly between the two 
limiting values of Q 2 . The rapid fluctuation of Ti then follows by subtraction of the BC 
value of T 2 . The importance of a' LT for a complete understanding of the spin structure has 
also been stressed by Li and Li ]18|. 



The GDH sum rule is obtained by integrating the quantity (o"i/2 — 0-3/2)/^ over the 
photon energy, where <t 3 / 2 and a 1 / 2 are the cross sections for virtual photon and target 
spins in parallel and antiparallel, respectively. Several authors have pointed out that serious 
discrepancies exist between the sum rule and theoretical analysis based on the existing 
unpolarized data. The small momentum evolution of the extended GDH sum rule was 
also investigated in the framework of heavy baryon ChPT Hl9fl . The authors predicted a 



positive slope of the GDH integral for Q 2 = 0, while the phenomeno logical analysis of 
Burkert et al. fl4| indicated a negative slope. However, the integrands of all spin-related 



sum rules are fluctuating functions of energy with both positive and negative contributions. 



Therefore, even seemingly small uncertainties may give rise to large differences (see Ref. [20 
and references quoted therein). 

The first direct experimental data have been recently taken at MAMI in the energy 
range 200 MeV< v < 800 MeV, and data at the higher energies are expected from ELSA 
within short. Concerning electroproduction in the resonance region, the A(1232) is being 
investigated at ELSA, Jefferson Lab, MAMI and MIT/Bates |2~2]] . The high-duty factor of 
these modern electron accelerators makes it possible to measure not only the (total) helicity 
cross sections but also the individual decay channels, thus providing a much more stringent 
test of model descriptions. The spin structure of the nucleon will be studied at the Jefferson 
Lab by a series of experiments ||23|| , which are expected to cover the entire resonance region 



and momentum transfers up to Q 2 ~ 2 (GeV/c) 2 . The outcome of these experiments should 
answer the open questions about the helicity structure at low and intermediate energies, the 
validity of the GDH sum rule, the size and role of the transverse-longitudinal contribution in 
connection with the BC sum rule and the generalized GDH integral, the evolution of these 
integrals with momentum transfer, and the contribution of individual decay channels to the 
spin structure functions. 

In the following we review the formalism of spin structure functions and sum rules in 
sect. 2. We then present the predictions of the UIM in sect. 3, and close by a brief summary 
of our findings in sect. 4. 



II. FORMALISM 



We consider the scattering of polarized electrons off polarized target nucleons. The 
lab energies of the electrons in the initial and final states are denoted with E and E', 
respectively. The incoming electrons carry the (longitudinal) polarization h — ±1, and the 
two relevant polarization components of the target are P z (parallel to the lab momentum 
k of the virtual photon) and P x (perpendicular to k in the scattering plane of the electron 
and in the half-plane of the outgoing electron). The differential cross section for exclusive 
electroproduction can then be expressed in terms of four "virtual photoabsorption cross 
sections" a^v, Q 2 ) by || 

da 



dtt dE 



7 = rV(z,,Q 2 ), (1) 



a = a T + ea L + hP x \jle{\ - e) a' LT + hP z y/l- e 2 cr' TT , (2) 

with 



a E' K 1 



(3) 



2tt 2 E Q 2 1 - e 

the flux of the virtual photon field and e its transverse polarization, v = E — E' the virtual 
photon energy in the lab frame and Q 2 > describing the square of the virtual photon 
four-momentum. In accordance with our previous notation || we shall define the flux with 
the "photon equivalent energy" K = fc 7 = (W 2 — m 2 )/2m, where W is the total cm energiy 
and m the mass of the target nucleon. 

The quantities ax and be expressed in terms of the total cross sections cr 3 / 2 

and 0i/2, corresponding to excitation of hadronic states with spin projections 3/2 and 1/2, 
respectively, 

= + 01/2) , O'rprp = —(03/2 01/2) • (4) 

The virtual photoabsorption cross sections in Eq. (Q) are related to the quark structure 
functions, 



K L v m 
^ = -^^ + »)' 

a TT = -^ (^1-7 <? 2 J , (5) 

with 7 = Q/v, and the quark structure functions Fx, F 2 , g\, and g 2 depending on v and Q 2 . 
In comparing with the standard nomenclature of DIS || we note that a' LT = —ctlt(DIS) 
and a' TT = —arT^DIS). 

We generalize the GDH sum rule by introducing the Q 2 -dependent integral 

2m 2 r x o I — \k 2 k for Q 2 — > 



g2 7o * ' + 0(Q- 4 ) for Q 2 - oo 



where x = Q 2 jlmv is the Bjorken scaling variable and xq = Q 2 / (2mm n + m 2 + Q 2 ) refers 
to the inelastic threshold of one-pion production. In the scaling regime the structure func- 
tions should depend on x only, and Tx — J gx(x)dx = const. We note that several other 
generalizations of the GDH integral have been proposed by introducing admixtures of the 
second spin structure functions g 2 and/or factors depending on (different) definitions of the 
photon flux. While all generalizations have the same limits as in Eq. (^j), they may differ 
substantially at intermediate Q 2 |25] . 



For the second spin structure function the Burkhardt-Cottingham (BC) sum rule asserts 
that the integral over g 2 vanishes if integrated over both elastic and inelastic contribu- 



tions [17]. As a consequence one finds 



2m 2 r 1 G M (Q 2 )-G E (Q 2 ) 2 

HQ ) = \ 92{x, Q ) dx = - — — — — — G M {Q ) , (7) 



Q 2 Jo * zv ' 4 1 + Q 2 /4m 2 

i.e. the inelastic contribution for < x < Xq equals the negative value of the elastic 
contribution given by the rhs of Eq. (|?p, which is parametrized by the magnetic and electric 
Sachs form factors Gm and Ge, respectively. These form factors are normalized by the 
nucleon's magnetic moment /ijv = kn + &n and its charge in the limit of Q 2 — > 0, 
Gm(0) = Hn an d Ge(0) = &n • The BC sum rule has the limiting cases 
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for Q 2 — > 



w) - < * . (8) 

C(g- 10 ) for Q 2 -> oo 



We have calculated the integrals ii and 7 2 in terms of the virtual photon cross sections, 

/ 2 g 2 = 5-2- / — - -a a 3/2 - a 1/2 --a' LT ) — , 10 

8ir 2 a Juo 1 + 7 2 \ 7 / f 

where i/q = m ?r + (^ 2 + Q 2 )/2m is the threshold lab energy of one-pion production. 

Since ^o' LT = 0(Q 2 ), the longitudinal-transverse term does not contribute to the integral 
Ii in the real photon limit. However, the ratio <j' lt /^ remains constant in that limit and 
hence contributes to I 2 . As a result we find 

J 2 (0) = -k 2 n + ^e N K N , (11) 

with the two terms on the rhs corresponding to the contributions of (73/2 — <J\i% and &' LT , 
respectively. Equation (|TT| ) provides an interesting model-independent constraint for the 
two spin-dependent virtual photon cross sections in the real photon limit. 

Finally, in order to relate our calculations with the experimental data, we express the 
virtual photon asymmetries A\ and A 2 defined in Ref. Q by the virtual photon cross sections, 

A 1 = ai/2 ~ ff3/2 =-^, A 2 = -^Z. (12) 
o"i/2 + 0"3/2 cr T cr T 

From the measured values for the asymmetries A\ and A 2 and from the total transverse cross 

section ax we then construct the quark structure functions gi and g 2 by use of Eqs. (f|), @ 

and (IE 



9l = -^(A l + 1 A 2 ), 92 = -^ (-A 1 + -A 2 ) . (13) 
1 + T 1 + T V 7 



In practice it is difficult to measure the two (unpolarized) cross sections <tt and gl separately 
by means of a Rosenbluth plot. The structure function F\ is then extracted from the total 
(unpolarized) cross section, a tot = o T + €Gl, by use of an "educated guess" for the ratio 



R = (JlI&t- Moreover, if only the longitudinal target polarization can be measured, only 
the combination A\ + 1]A 2 can be determined, where rj = eQ/(E — eE'). In this case the 
structure function g\ has to be extracted by use of additional assumptions about the second 
spin structure function g 2 . 



III. RESULTS AND DISCUSSIONS 



A. Total cross sections 



The recently developed Unitary Isobar Model (UIM) || describes the single-pion electro- 
production channel quite adequately. However, for energies W > 1.3 GeV the contributions 
from other channels, in particular from multipion production, become increasingly impor- 
tant. At present several theoretical models are available for the two-pion photoproduction 



channels [26— 28| , but little information exists on multipion electroproduction. In this paper 
we extract the necessary information from the available data for the total cross sections 
of real and virtual photon absorption on the proton ]2P,50|, assuming that the single-pion 
channel is well determined by the UIM. In addition we account for the well investigated eta 
production channel by using the parametrization for the corresponding total cross section 
(W, 0) at the photon point as suggested in Ref. ||10|| . We extend this parametrization for 



the case of virtual photons using the helicity amplitude A\^ 1 2 1 for the 5*11(1535) resonance 



obtained within the UIM, i.e. 



dSu), 



(14) 



The transverse part of the total cross section for multipion production is now obtained 



by 



4"° W Q 2 ) = F 2 (Q 2 ) a¥ xp, (W, 0) - a { T L7T, (W, 0) - a™(W, 0) 



(exp) 



(It), 



r (l), 



(15) 



with a monopole form factor F(Q 2 ) = 1/(1 + a nn Q 2 ) describing the Q 2 dependence. A 
value of a nn = 1.05 (c/GeV) 2 provides the best fit to the total inelastic cross section up to 



the third resonance region and for Q 2 < 1 (GeV/c) 2 . In Fig. 1 we present our results for the 
total cross sections at Q 2 = and 0.5 (GeV/c) 2 . 

We have not yet included the contributions of the eta and multipion channels to the 
longitudinal part of the cross section. Since the ratio R is small in the resonance region, 
R = ai/ 'ax ~ 0.06 PD|,|BJ, these contributions should not be of major importance anyway. 



On the other hand, we find a substantial longitudinal strength for single-pion production 
at the higher energies, R {ln) = 4^V4^ ) > 0.5 at W > 1.8 GeV and Q 2 > 0.3 [GeVjcf. 
Hence we are led to the conclusion that the (small) longitudinal contribution is dominated 
by single-pion production. 

B. Asymmetries and structure functions 

One of the basic ingredients for the asymmetries A\ i and structure functions gi^ is the 
difference of the two helicity cross sections, 

Aa = ct 3/2 - en/2 = Ao- (lw) + Aa M + A<r (n7r) . (16) 

The single-pion contribution to this difference has been calculated within the UIM ap- 
proach. Since the dominant mechanism of eta electroproduction is due to the excitation of 
the Sii(1535) resonance, the contribution of the eta channel to Eq. (|16|) may be approxi- 
mated by Aa^ « — 2<jrp . The description of the multipion channels is the most delicate 



part in our calculation of the sum rules, because the mentioned theoretical models [f2q -|28 
for two-pion photoproduction were tuned to the existing data for total cross sections. The 
simplest estimate for A<T ( - n7r - ) is based on the assumption that the two-pion contribution is 
generated by resonances, and that its helicity structure follows the known behavior of the 
one-pion contribution Ac^ 17 ^ pl|j3^|, i.e. Acx^ n7r - ) /^ cr '' l7r ^ = const. In the present paper we 



first account for the one-pion and eta channels explicitly, and then use the simple prescription 

(nTr) 

Act'™' = ^— (AaW + Aa^) , (17) 
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with the factor a^ 1 in the numerator providing the correct threshold behaviour for the 
two-pion contribution. 

In the upper part of Fig. 2, our results are compared to the asymmetry Ai+r]A 2 measured 
at SL AC [[J . As may be seen, the prescription of Eq. (|TT| ) is in reasonable agreement with 
the data up to W = 2 GeV. We also note that the contribution of the r\ channel (dotted 
curves) leads to a substantial increase of the asymmetry over a wide energy region. 

Having fixed all the ingredients we can now calculate the structure functions and gener- 
alized GDH integrals. Our results for the structure function gi are presented in the lower 
part of Fig. 2. Up to a value of W 2 = 2 GeV 2 (corresponding to x = 0.31 and x = 0.52 at 
Q 2 =0.5 and 1.2 (GeV/c) 2 , respectively), the main contribution to g± is due to single-pion 
production. The negative structure above threshold is related to excitation of the A (1232) 
resonance. In the second and third resonance regions the contributions from the eta and 
multipion channels become increasingly important. Note that with increasing value of Q 2 
the relative importance of the eta channel increases. This phenomenon is connected with 
the well known peculiarity that the helicity amplitude of the Sii(1535) resonance decreases 
much slower with Q 2 than for any other resonance f34||9|1. 

C. Integrals I\ and I 2 

In Fig. 3 we give our predictions for the integrals I\{Q 2 ) and ^(Q 2 ) i n the resonance 
region, i.e. integrated up to W max = 2 GeV. The striking feature of the integral l\ is 
the evolution from large negative to small positive values in order to interpolate from the 
GDH prediction for real photons to the data obtained from DIS. As can be seen from the 
top of Fig. 3, our model is able to generate the dramatic change in the helicity structure 
quite well. While this effect is basically due to the single-pion component predicted by the 
UIM, the eta and multipion channels are quite essential to shift the zero-crossing of l\ from 
Q 2 = 0.75 (GeV/c) 2 to 0.52 (GeV/cf and 0.45 (GeV/c) 2 , respectively. This improves the 
agreement with the SLAC data || . However, some differences remain. Due to a lack of data 
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in the A region (see lower part of Fig. 2), the SLAC data are likely to underestimate the A 
contribution. Since this contribution is negative, a numerical integration of the data points 
could overestimate the I\ integral or the corresponding first moment Y\. A few more data 
points in the A region would be very useful in order to clarify the situation. 

Concerning the integral I2, our full results are in good agreement with the predictions 
of the BC sum rule (see Fig. 3, lower part). The remaining differences are of the order 
of 10 % and should be attributed to contributions beyond W max = 2 GeV and the scarce 
experimental data for a' LT . With regard to the latter problem we note that the BC sum 
rule, Eqs. (|7|) and (|H]), contains very sizeable contributions from the longitudinal-transverse 
cross section. In Table 1 we list the contributions of the different ingredients of our model 
to the integrals Ji and I2 at the real photon point. 

The convergence of the sum rules can not be given for granted. In fact Ioffe et al. |33| 



have argued that the BC sum rule is valid only in the scaling region, while it is violated by 
higher twist terms at low Q 2 . Therefore the good agreement of our model with the BC sum 
rule could be accidental and due to a particular model prediction for the essentially unknown 
longitudinal-transverse interference term. As can be seen from Fig. 3, the contribution of 
ctlt' is quite substantial for I2 even at the real photon point due to the factor 7 _1 in Eq. 
10. This contribution, however, is constrained by the positivity relation \<j' lt \ < ^Jol&t = 
\J~RaT- The dash-dotted line shows the integral for the upper limit of this inequality and 
a similar effect would occur for the lower limit. This surprisingly large contribution can be 
understood in terms of multipoles. In a realistic description of the integrated cross section 
<jlt' the large M\ + multipole can only interfere with the small L\ + multipole. The upper 
and lower limits of the positivity relation overestimate the structure function considerably 
due to an unphysical "interference" between s- and p-waves. 

In Fig. 4 we compare our result for the function T\, defined in Eq. (6), with the 
calculation of Ioffe []15| for the resonance contribution. This representation demonstrates 
more clearly the importance of the higher channels, e.g. 77 and 2ir. 

Fig. 5 shows our predictions for Aa = cr 3 / 2 — cr 1 / 2 at Q 2 =0, 0.5 and 1.2 (GeV/c) 2 . 
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In particular we note the following peculiarity of the Q 2 dependence of the resonance and 
multipion contributions. With increasing value of Q 2 the positive contribution of the A 
resonance decreases rapidly due to the strong NA transition form factor. In the second and 
third resonance regions, however, Act changes from positive to negative values of Q 2 . The 
reasons for this zero-crossing are (I) the helicity structure of the -D13 and F 15 resonances, 
which changes from cr 3 / 2 to a\/2 dominance f35[, and (II) the surprisingly slow fall-off of the 
S\\ transition form factor |34||| , which increases the negative contribution of this resonances 
relative to the others. 

Both effects conspire to cancel the A contribution already at Q 2 w 0.5(GeV/c) 2 and lead 
to the observed positive value of I\ at larger momentum transfer. The same cancellations 
are responsible for the rapid decrease of I2 with Q 2 . However, it should be kept in mind 
that this integral also receives strong contributions from u' LT . In view of our bad knowledge 
of the longitudinal-transverse response we find it actually quite surprising that the UIM 
"knows" about the BC sum rule and the predicted rapid fall-off with Q 2 . 



IV. CONCLUSION 

We have studied the nucleon spin structure functions and the generalized GDH and BC 
integrals for small and moderate momentum transfer within the framework of the Unitary 
Isobar Model 0. The contributions from eta and multipion channels in the second and 
third resonance regions have been modeled by accounting for the total photoproduction 
cross section. 

We find that the role of the eta channel becomes more and more important with increas- 
ing values of Q 2 . At the real photon point, this channel tends to cancel the contribution of 
multipion production. However, at finite Q 2 the contributions of both channels add coher- 
ently. As a result the zero-crossing of the integral I\ shifts from Q 2 = 0.75 (GeV/c) 2 (for the 
one-pion channel only) to Q 2 = 0.45 (GeV/c) 2 (if eta and multipion channels are included). 

Our analysis indicates that both the measured asymmetries and the theoretical models 
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have to be quite accurate in order to determine the integrals I\ and I2 in the resonance 
region. As far as theory is concerned, a more quantitative description of the multipion 
channels is absolutely necessary. On the side of the experiment, a measurement of all four 
virtual photon cross sections is needed to determine the quark spin structure functions 
in a model-independent way. Moreover, the oscillatory behaviour of the integrand in the 
resonance region will require small error bars and a relatively small energy binning in order 
to give precise values for the GDH and BC integrals. The proposed experiments at Jefferson 



Lab |23[ carry the promise to obtain such data in the resonance region and we are eagerly 
waiting for the outcome of these investigations. 
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TABLES 



h,2 


Born+A 


Pu,D 13 ,... 


V 


multipion 


total 


sum rule 


h 


-0.565 


-0.152 


0.059 


-0.088 


-0.746 


-0.804 


h 


1.246 


0.063 


-0.059 


0.088 


1.338 


1.252 



TABLE I. Contributions of the different channels to the integrals I\ and I2 at the real photon 
point, Q 2 = 0. The column "Pn,Di3, ..." shows the contributions of all the resonances above the 
A that are included in Ref. |J. Note that the column "total" lists the values for integration up to 
Wmax = 2 GeV only, while "sum rule" is the prediciton for the full integral. 
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FIGURES 
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FIG. 1. Total cross sections for photoabsorption and inelastic electron scattering on the proton 
for e = 0.9 and Q 2 = 0.5 (GeV/ c) 2 . Dashed, dotted and dash-dotted curves: contributions of 
single-pion, eta and multipion channels, respectively; solid curves: final result. Experimental data 
for the total cross sections from Refs. p(| (x) and |29| (o), for the two-pion production channels 
from Ref. |29) (A). 
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FIG. 2. The asymmetry A\+rjA2 (top) and the spin structure function g\ (bottom) as function 
of the Bjorken scaling variable x at Q 2 = 0.5 and 1.2 (GeV/c) 2 . Dashed, dotted and solid curves: 
calculations obtained with In, lir + rj, and lTr + rj + nir contributions, respectively. Data from Refs. 
[||] (•) and ]36| (o). The error bars give the statistical errors, the systematical errors are estimated 
to be of equal size. 
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FIG. 3. The integrals I\ and I2 defined by Eqs. (10) and (11) as functions of Q 2 in the 
resonance region, integrated up to W max = 2GeV. Upper figure: Notation as in Fig. 2 and data 
from Ref. ||. Lower figure: The full and dashed lines are our predictions with and without a' LT 
(see Eq. (10)), the dash-dotted line is obtained for a' LT = ^J<jl&t (see text), and the dotted line 
is the sum rule prediction of Ref. [1 71]. All calculations for I2 include lir + r/ + nir contributions. 
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FIG. 4. The integral T± for the proton defined by Eq. (6) as function of Q 2 . The solid and 
dashed lines show the integral up to W max = 2GeV for our full calculation and our calculation 
without the rj and 2tt contributions, respectively. The dash-dotted line is obtained from Ref. [15] 
and contains only the resonance component. The data are from Ref. ||. 
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